Breakdown of time-dependent mean-field theory for a 
one-dimensional condensate of impenetrable bosons 
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We show that the time-dependent nonlinear Schrodinger 
equation of mean-field theory has limited utility for a one- 
dimensional condensate of impenetrable bosons. Mean-field 
theory with its associated order parameter predicts interfer- 
ence between split condensates that are recombined, whereas 
an exact many-body treatment shows minimal interference. 
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Mean-field theory (MFT) has proven remarkably suc- 
cessful at predicting both the static and dynamic be- 
havior of Bose-Einstein condensates (BECs) in weakly- 
interacting atomic vapors [jjj, including the ground state 
properties |2[q], the spectrum of collective excitations 
{§[i|, four- wave mixing HQ, matter- wave solitons |§J§, 
and interference between BECs jl(i|-|l"5)|. The basic no- 
tion underlying MFT is of a macroscopic wave function 
jl4],[t5]|, or order parameter, which defines the spatial 
mode into which a significant fraction of the atoms con- 
dense below the critical temperature. The macroscopic 
wave function typically obeys a nonlinear Schrodinger 
equation (NLSE), the Gross-Pitaevskii equation, and is 
most suitable for dilute Bose gases. But the success 
of MFT is not assured in all cases. For example, in 
one-dimension (ID) a spatially homogeneous ideal gas 
in its many-body ground state exhibits complete BEC 
into the lowest single-particle state, but no BEC at any 
nonzero temperature. Furthermore, previous exact anal- 
ysis [|l6|,[nj of a gas of impenetrable or hard core bosons 
in ID by one of us (MG) and its extension by Lenard 
p8[ and Vaidya and Tracey jl9| have shown that in the 
many-body ground state the occupation of the lowest 
single-particle state is of order ^/~N where N is the total 
number of atoms, in contrast to N for usual BEC. Nev- 
ertheless, since N 3> 1 and the momentum distribution 
has a sharp peak in the neighborhood of zero momentum 
Ipjl^ot , this system shows some coherence effects. More- 
over, Olshanii [^0| has shown theoretically that such a ID 
gas of impenetrable bosons can be realized at sufficiently 
low temperatures and densities in thin atom waveguides, 
and there is also considerable experimental effort towards 
realizing atom waveguides for producing ID atom clouds 
pl| , p2t , a new and fertile system for studying the physics 
of condensed systems. 

Our goal in this Letter is to show that mean-field the- 
ory breaks down in the analysis of the dynamics of ID 



atom clouds, in that it predicts interference effects that 
are absent in the exact theory. In particular, we have re- 
cently analyzed many-body solitons in a ID gas of bosons 
using the time-dependent generalization of the Fermi- 
Bose mapping p3| , which maps a ID gas of hard core 
bosons to a gas of free fermions jl6],|4j , and Kolomeisky 
et al. [P5[ have proposed a corresponding NLSE with a 
quartic nonlincarity to extend the usual mean-field the- 
ory for ID atom clouds. For a harmonic trap the ground- 
state density profiles from their theory show excellent 
agreement with the exact many-body results (see Fig. 1 
of their paper). The key question, then, is whether this 
extended NLSE can be used in all circumstances. To ad- 
dress this issue we examine the problem of a ID atomic 
cloud in the ground state of a harmonic trap that is 
split by a blue-detuned laser and recombined, both using 
an exact many-body treatment based on the Fermi-Bosc 
mapping and the approximate NLSE: the NLSE predicts 
interference whereas the exact analysis does not. 

The fundamental model we consider is a ID gas of N 
hard core bosons at zero temperature described by the 
Schrodinger Hamiltonian 
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where Xj is the one-dimensional position of the jth par- 
ticle and the many-body potential V is symmetric (in- 
variant) under permutations of the particles. The two- 
particle interaction potential is assumed to contain a hard 
core of ID diameter a. This is conveniently treated as a 
constraint on allowed wave functions ip{x\, • • • , xn', t): 

ip = Q if \xj-x k \<a , l<j<k<N, (2) 

rather than as an infinite contribution to V, which 
then consists of all other (finite) interactions and ex- 
ternal potentials. To construct time-dependent many- 
boson solutions of Eq. (|l|) we employ the Fermi- 
Bose mapping |l^,|2^,|24| and start from fermionic so- 
lutions ip F (xi, ■ ■ ■ , xn; t) of the time-dependent many- 
body Schrodinger equation (TDMBSE) Hvb = ihdip/dt 
which are antisymmetric under all particle pair exchanges 
Xj ^ ' X • hence all permutations. Next introduce a "unit 
antisymmetric function" 

A(xi, ■ ■ ■ ,x N ) = Yl s E n ( x k 

l<j<k<N 
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where sgn(a?) is the algebraic sign of the coordinate dif- 
ference X — X X j •} 1.6. j it is +1(-1) if x > Q(x < 0). For 
given antisymmetric ipp, define a bosonic wave function 

ijj B by 

ip B (xi, - ■ ■ ,x N ;t) = A(xx, - ■ ■ ,x N )ip F (xi, - ■ ■ ,x N ;t) (4) 

which defines the Fermi-Bosc mapping. Then tps sat- 
isfies the hard core constraint (2) if i/jp does, is totally 
symmetric (bosonic) under permutations, and obeys the 
same boundary conditions In the Olshanii limit 

pOf (low density, very thin atom waveguide, large scat- 
tering length) the dynamics reduces to that of the impen- 
etrable point Bose gas, the a — > limit of Eq. (|^). Then 
under the assumption that the many-body potential V of 
Eq. (|I|) is a sum of one-body external potentials V(xj, t), 
the solution of the fermion TDMBSE can be written as 
a determinantal wavefunction l|lq,E4| 
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V>f(xi, • • • , x N ; t) — C dct <t>i(xj,t), 
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where 
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It then follows that ipF satisfies the TDMBSE, and it 
satisfies the impenetrability constraint (vanishing when 
any Xj = xg) trivially due to antisymmetry. Then by 
the mapping theorem ipB of Eq.(|^) satisfies the same 
TDMBSE. For each i = 1, . . . , N the initial normalized 
wave functions 4>i{x, 0) represent different orbitals of the 
potential at t = 0. For example, for a harmonic trapping 
potential (f>i(x, 0) are given by the well-known Hcrmitc- 
Gaussian orbitals, with one particle occupying each of 
the N lowest levels for the iV-particle ground state p5| . 

We have performed numerical simulations based on 
Eqs. (JsJ) , using the split-step beam propagation method 
p6[ , for an initial harmonic trap containing N particles 
in their ground state which is subjected to a centered 
Gaussian repulsive potential, e.g., a blue-detuned laser 
field which serves to split the initial state. After some 
time tp t both the harmonic trap and repulsive poten- 
tial are turned off and the two split components allowed 
to recombine: this is an interference experiment of the 
cool, cut, interfere variety previously discussed jllj]. The 
potential in Eq. (0) is taken of the specific form 



V(x, t) = -muo^x + Vb sin( 
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and V = for t > t pot . Figure shows an illustrative 
example for N — 10 with a repulsive potential of height 
Vb = 20huj, and width w = 3xo, with xq — y/h/2muj 
the ground state harmonic oscillator width, and Lot pot = 
3. This figure shows a gray scale plot of the TV-particle 



density pn (x, t) = \4>i{ x ^)\ 2 |l6| , p3| , p4| as functions 



of Lot (horizonatal axis) and position x/xq (vertical axis), 
with white being the highest density. 




FIG. 1. Exact many-body theory simulation of the cool, 
cut, interfere scenario. The figure shows a gray-scale plot 
of the particle density pjf(x,i) as a function of uit (horizon- 
tal axis) and position x/xq (vertical axis), with white being 
the highest density, for = 10, Vb = 20hu, w — 3xq, and 



pot 



= 3. 



The potential height was chosen such that Vb > A 4 = 
Tiuj(N — 1/2), where p, is the chemical potential of the 
TV-particle oscillator ground state GjJ, noting that the 
top of the iV-particle Fermi sea is at n — N — 1. As ex- 
pected, as the repulsive potential turns on it splits the ini- 
tial ground state into two separated components. Upon 
release at t = t pot the two components expand and subse- 
quently recombine. What is noteworthy is that although 
there is some modulation upon recombination there are 
no strong modulations indicative of the interference pro- 
visionally expected for bosons: this was a generic find- 
ing from our simulations irrespective of the time scale 
on which the repulsive potential was turned on |0,tt2[. 
In contrast, the density |^>i(a;,£)| 2 for each individual or- 
bital i = 1, . . . , N can show large inteference fringes, but 
with a different period in each case. Thus, the minimal 
inteference seen in Fig. 1 is a result of washing out of 
the individual interferences by averaging over N orbitals. 
Thus, the remnant of any interference fringes decreases 
with increasing N and vanishes in the thermodynamic 
limit. 

Physically, it makes sense that the interference fringes 
are all but absent since the Fermi-Bose mapping shows 
that in this ID limit the system of bosons acts effectively 
like a system of free fermions insofar as effects expressible 
only in terms of |?As| 2 are concerned, so interference is 
not expected ^7j. The lack of interference is therefore 
a signature of the Fermi-Bose duality that occurs in ID 
systems of impenetrable particles fl6|| . 

We next turn to the mean-field description proposed 
by Kolomeisky et al. [ p5| for low-dimensional systems. 
In particular, they introduce an order parameter $(x,t), 
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normalized to the number of particles N, for such sys- 
tems, though they do not discuss what this represents 
physically. Using energy functional arguments they de- 
duce the following NLSE with quartic nonlinearity for a 
ID system of impenetrable bosons: 



lh ~dt = 
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with V(x, t) the same as in Eq. (||). Our goal now is to 
compare the predictions of this NLSE for the same cool, 
cut, interfere simulation in Figure [l], with the initial con- 
dition <j>(x,0) = yj pn(x, 0) corresponding to the exact 
many-body solution, all other parameters being equal. 
Figure || shows the corresponding gray-scale plot of the 
density p(x,t) = \Q(x,t)\ 2 , and two features are appar- 
ent: First, during the splitting phase when the repulsive 
potential is on there is very good overall agreement be- 
tween the exact many-body theory and the NLSE predic- 
tion. Second, when the split components are released and 
recombine they produce pronounced interference fringes, 
in contrast to the exact theory. Indeed, this interference 
in the MFT is to be expected on the basis of previous 
theoretical work [jl0| , even though a quartic (rather than 
quadratic) nonlinearity is employed here. Thus the MFT 
cannot accurately capture the time-dependent dynamics 
in all situations. 




FIG. 2. Mean-field theory simulation of the cool, cut, in- 
terfere scenario. The figure shows a gray-scale plot of the 
particle density p(x,t) — j$(a;,t)| 2 as a function of u>t (hori- 
zonatal axis) and position x/xq (vertical axis). Parameters 
are the same as Fig. 1. 

In conclusion, the numerical simulations presented 
show that MFT seriously overestimates the magnitude 
of interference phenomena in a ID gas of impenetrable 
bosons. However, Kolomeisky et al. have demonstrated 
that their NLSE produces excellent results for the ground 
state density profile of ID atomic clouds, and our sim- 
ulations show that during the splitting process there is 
excellent agreement between the exact theory and the 



MFT. During the splitting the density profile is found to 
closely follow the ground state of the local applied poten- 
tial, harmonic trap plus repulsive potential, meaning that 
the phase of the order paramter is not a key element dur- 
ing this time interval. In this case the order parameter 
for the NLSE is physically Q(x,t) = \/p(x, t), to within 
an arbitary overall phase factor. In contrast, during the 
recombination phase the relative phase of the order pa- 
rameter between the two split components is the key 
element allowing for the appearence of the interference 
fringes in the MFT. Since the exact many-body theory 
does not show such fringes this means that the MFT ap- 
proach endows the order parameter with phase informa- 
tion that is beyond the real degree of coherence present 
in the system. Thus, as long as this coherence does not 
manifest itself in the form of interference, the MFT is 
of some utility, but not if the predictions are sensitive 
the the order parameter phase. This view is consistent 
with the lack of off-diagonal-long-range-order in the ID 
gas of impenetrable bosons jl8|,[l9) , which is normally re- 
quired for the introduction of an order parameter Jl4|,|l5| . 
In future work we hope to investigate the transition of 
impenetrable bosons from the ID regime, with no inter- 
ference effects, to the 3D regime, which is well described 
by MFT and shows interference jt(^. In particular, at 
ID atom waveguide junctions such as atom beam split- 
ters or recombiners, the system will be sensitive to the 
fact that the ID waveguides are actually embedded in a 
higher-dimensional space, in which case atom interferom- 
eters employing ID atom waveguides may still be capable 
of producing interference fringes. Another very interest- 
ing question is whether or not the usual GP NLSE or its 
2D modification proposed by Kolomeisky et al. j2jj] gives 
reliable predictions of interference between BECs in the 
2D regime. This has potentially important experimen- 
tal consequences in view of progress in construction of 
2D atom waveguides and their potential applicability to 
atom interferometry p8fl . 
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